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In the orbitally degenerate ( J = 5/2 ) Periodic Anderson Model, the mag- 
netic susceptibility is composed of both the Pauli term and the Van Vleck 



term, as is well known. The former is strongly enhanced by the strong cor- 



relation between /-electrons. But, for the latter, the influence of the strong 
correlation has been obscure for years. In this paper we give the solution of 

•l-H . 

^ ■ the longstanding problem. With the aid of the d = oo approximation, we 

study this problem on the basis of the Fermi liquid theory with degenerate 
orbitals, taking account of all the vertex corrections in a consistent way. As 
a result, we obtain the simple expression for the magnetic susceptibility, and 
show unambiguously that the Van Vleck term is also highly enhanced in the 
strong correlation regime. This fact explains naturally the enhanced magnetic 
susceptibility observed in many insulating systems ( i.e., Kondo insulator ). 
Moreover, we show that the Wilson ratio takes a value around 1 in the metallic 
system, in good agreement with experiments. 
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I. INTRODUCTION 



In this paper, we investigate the uniform magnetic susceptibility, Xi f° r the orbitally 
degenerate ( J = 5/2 ) Periodic Anderson Model ( J = 5/2 PAM ). [0,0 This model is 
much closer to real heavy Fermion systems because it includes both the /-electron orbital 
degeneracy and the l-s coupling. Based on this model, we can explain several characteristic 
behaviors observed experimentally, which are unable to be derived from the SU(N)-PAM. ( 
SU(N)-PAM is familiar for theorists but less realistic model than J = 5/2 PAM in that both 
the /-electron bands and the conduction electron bands possess the same N-fold degeneracy. 

) 

One of the characteristic behaviors of the model is the presence of the anomalous Hall 
effect (AHE). Reference |§ clearly explains that the anomalous Hall coefficients of heavy 
Fermion systems are proportional to the square of the resistivity at low temperatures, and 
points out that its mechanism is similar to that of AHE on the ferromagnetic metals proposed 
by Karplus and Luttinger. j|,|5| 

The other characteristic property is the presence of the Van Vleck magnetic susceptibility. 
H|7| As is well known, the magnetic susceptibility of this model contains the Van Vleck 
susceptibility xv in addition to the usual Pauli susceptibility xp-, so the total magnetic 
susceptibility is given by x = Xp + Xv- Xp comes from the intra-band contribution, which 
vanishes in the absence of the Fermi surface, and xv comes from the inter-band contribution. 
Of course, x is one °f the most important and characteristic physical quantities in the 
heavy Fermion systems. But, unfortunately, the property of the Van Vleck susceptibility 
Xv remains quite obscure in the presence of the large Coulomb interaction. 

Furthermore, its general properties should be studied for understanding the origin of the 
strongly enhanced magnetic susceptibility for so-called Kondo insulators, f| where xp = 
and only the Van Vleck susceptibility remains. We also have to understand the property of 
X in order to interpret the Knight shift data for super-conductive heavy Fermion systems. 
1 

Historically, this problem was first studied in ref. Q. The authors paid attention only 
to xp-, an d claimed that xv could be neglected. On the contorary, mean-field like treat- 
ments give the enhanced xv ■ [l^UHl But, such approximations seem to be insufficient to 
settle this problem because they assume the frequency independent renormalization factor, 
which is a less-realistic assumption. The settlement of this problem will bring us the useful 
understanding for the electronic properties for the orbitally degenerate model. 

The existence of xv i n J — 5/2 PAM is closely related with the fact that the magne- 
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tization operator M is not conserved in this model, i.e., [H,M] ^ 0. [12] Whenever M is 
conserved, xv vanishes identically and the susceptibility has been obtained without ambi- 
guity in the framework of the Fermi liquid theory ( FLT ). For instance, the susceptibility 
for SU(2)-PAM or J = 5/2 single-site Anderson model is made of xp alone, and obtained by 
the Luttinger's method. ||13|-4T5"f On the other hand, in order to analyze the magnetic suscep- 
tibility of our model, we have to find another useful method. ( By the RPA-approximation, 
which is an insufficient analysis for our problem, the enhancement for xv is l ess than factor 
of two even in the strong correlation regime where xp ~ 00 • [0 ) 

In this paper, we analyze the properties of x an d Xv by investigating all the vertex 
correction for them by use of a kind of the Bethe-Salpeter equation. Then we employ the 
d = oo approximation, JIT] which is considerably effective for our problem on J = 5/2 PAM. 
rj| We also discuss the Wilson ratio R in the strong correlation regime. The definition of 



R is 

where R° is the unrenormalized Wilson ratio. 7 ( 7 ) is the ( unrenormalized ) T-linear coef- 
ficient of the specific heat, x ( X° ) is the ( unrenormalized ) uniform magnetic susceptibility. 
R° is given by 

do _ 3x° 

2sVUV°(o)' 1 } 

where p°(0), g and \xb are the density of states for /-electrons, the Lande's ^-factor, and 
the Bohr magneton respectively, and J^ s = J (J + 1) in the spherical system, i.e., without 
any electronic crystal field (ECF). For instance, for an isolated Kondo atom R is universally 
given by R = (2J + 1)/2J and R° = 1. 

The composition of this paper is as follows : in §2, we explain main properties of J = 5/2 
PAM and its Green's functions. In §3, we express the magnetic susceptibility in terms of the 
dynamical susceptibility on the basis of the multicomponent FLT. In the insulating systems, 
the susceptibility is given by only the Van Vleck susceptibility defined in this section. In §4, 
we calculate x an d Xv explicitly in case U = 0. The obtained results are the same as those 
obtained previously. In §5, we investigate x an d Xv(p) m d = 00 in the strong correlation 
regimes, by taking all the vertex corrections into account. Then, in the case of no ECF, we 
find that both x and xv are proportional to the same enhancement factor. This implies that 
the Van Vleck susceptibility xv is enhanced, which is the main result of this paper. In §6, we 
also discuss the Wilson ratio in our system and find out that R ~ 1. We also point out that 



3 



the results for d = oo systems are qualitatively correct for the three dimensional system. In 
§7, the obtained results are summarised and the future problems are pointed out. 

In Appendix A, we study the magnetic susceptibility for SU(N)-PAM, and show that 
the Van Vleck susceptibility is small in this model, which arises from the difference between 
the magnetic momentum of /-electrons and that of conduction electrons. In Appendix 
D, we confirm our conclusions numerically, by use of the ( self-consistent ) second order 
perturbation treatment ( ( SC- ) SOPT ). In Appendix E, we show in detail that given 
by (fT8|), vanishes identically. 



II. MODEL AND GREEN'S FUNCTION 



In this paper, we study J = 5/2 PAM, which is a realistic model for the Ce-compound 
heavy Fermion systems in that it represents the /-orbital degeneracy. In the absence of the 
magnetic field H, our Hamiltonian is given by 0|| 

H = H + H 1 (3) 

ko- kAf A/ko- 

H\ = — ^2 /k-qA//k'+qA/'/k'M'/kM, 
Z kk'qM^M' 

where cj^ is the creation operator of a conduction electron with momentum k and spin er, 
and fljM is the creation operator of an /-electron with momentum k and angular momentum 
M, which is the eigenvalue of J z = l z + s z . As J = 7/2 states are ignored in our model as 
is usually done, the /-orbital has six-fold degeneracy, i.e., M = {5/2, 3/2, • • • , —5/2}. VkMo- 
in (f|) is the mixing potential between the f-electrons and the conduction electrons, which is 
given by 

W = (4vr)V2 £ a&Y&iGk, Vk ) ■ V, (4) 

m 

aZ = -a{(7/2 - Ma)/7} 1/2 5 m , M ^ /2 for J = 5/2, 

where a^ a is the Clebsh-Goldan ( C-G ) coefficient and Y™ 3 (9k, <fh) is the spherical harmonic 
function. In the Coulomb interaction term @, we take account of the Pauli's principle, i.e., 
M ^ M'. 

In the presence of the magnetic field H along z-axis, both and in (|) are shifted 
by the Zeeman energies as 
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E f M = Ef + gf i B M-H, 

(5) 

Cka = Ck + ^B^" • H, 

where li-q is the Bohr magneton and a = 1 for up-spin and a = — 1 for down-spin, g is the 
Lande's ^-factor ( g = 6/7 for J = 5/2 ). Hereafter, we put \i& = 1. 

The general forms of the Green's functions for J = 5/2 PAM in the absence of the 
magnetic field are given by (3.11) ~ (3.13) of ref. |§. In the presence of the magnetic field 
H, the /-electron Green's functions can be expressed as 

GkMM'M = G ImM'( U ) + akAfaM G k CTCT '(^) a kAf'<x'M; ( 6 ) 

GImm>(<») = (wl Sk(w))^ # , (7) 

Q^kMo- = y^^kM'o-^kM'M^); (8) 
AP 

where /x denotes the Fermi energy of this system and G-ka- a r(u) is the conduction electron 
Green's function. Considering the Zeeman terms for only the /-electrons, we can derive the 
expression for the conduction electron Green's function in the presence of the magnetic field 

q ,s <W • (^ + V ~ e k) + (1 ~ 2(5 <T(T /)S' w (u;) 

+ - e k - S kn (u))(uj + li - e k - S m (u)) - S^^S^uj) ' 

where 

= I] V ±MaGl MMI {u)V£ MltT ,. (10) 
MM' 

Note that Sk CT ,- CT (u;) = in case H = 0. The pole of Gk CT (a;) represents the spectrum for the 
quasiparticle ( i.e., heavy electron ), which we denote as ££. The pole of G^{u) represents 
the local /-electron spectrum, , and never represent the quasiparticle spectrum. GrkM (<*■>) 
contains both poles. 

In Fig. 1, the electronic structure of J = 5/2 PAM in case U = is sketched in terms 
of the one-body picture. There remains four-fold degenerate local /-electron spectrum, 
which reflects the difference of the degeneracy between the conduction electrons and the 
/ electrons. ( They do not exist in SU(N)-PAM. ) In Fig. 1(a), the Fermi energy \x lies 
below the lower edge of the f-c hybridization gap, so the system is metallic, kp is the 
Fermi momentum and D ( — D ) is the upper ( lower ) edge of the conduction band. We 
can show that V 2 / D 2 is proportional to ^/-band/^c; where mj-b an d and m c represent the 
unperturbed mass of f-electrons and the conduction electrons, respectively. In usual heavy 
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fermion systems, V 2 / D 2 <C 1 is satisfied. [18 , 19[] On the other hand, in Fig. 1(b), the 
Fermi energy \x lies in the hybridization gap and below E*, so the system is insulating. 
A_ ( A + ) represents the lower ( higher ) edge of the hybridization gap measured from 
the Fermi energy. This situation is a prototype for the so-called Kondo insulator. [|5] More 
detailed electronic structure is studied beyond the one-body picture by use of the numerical 



perturbation calculation. \TE 

In the absence of the magnetic field, Gkaa'i^) is diagonal with respect to the spin a 
even if U ^ 0. || The density of states ( DOS ) of the /-electrons, p°(uj), and that of the 
conduction electrons, p c °(uj), are given by the retarded Green's functions as 

PcV) = -^£hnG°», (11) 

k 

PV) - -^X>«(") = {u + f_ Ef? ■ AM, (12) 



where N represents the total number of the /-sites. In case U ^ 0, these relations (|TT|) and 
(|12|) are also satisfied with E* replaced by E* + £(&>) for \u\ < T*, where T*(> 0) represents 
the characteristic energy within which the quasiparticles are well defined. In usual heavy 
Fermion systems, T* is much smaller than D. 

As is shown in ref. ||, G^mm 1 ^) oc e^ M '~ M ^ k even in case U ^ 0, which is an important 
relation throughout this paper. 



III. DEFINITION OF THE PAULI AND VAN-VLECK SUSCEPTIBILITY 

In this section, we consider both the Pauli susceptibility xp an d the Van Vleck suscep- 
tibility xvi an d derive their general expressions on the basis of the FLT for J = 5/2 PAM. 
Though our discussion in this section is restricted to J = 5/2 PAM ( in the presence of any 
ECF ) for simplicity, it is more general and is valid for many kinds of multicomponent Fermi 
liquid systems. Now, we consider the situation where the system is in a finite magnetic field 
H along 2-axis. The magnetization of this system, (M), is expressed by use of the Green's 
functions as 

( M > = A?E £ ^-e +i -°|E^MG k MMM + E^ kCT(7 M|. (13) 

The second term of the above equation is the contribution from the spin of the conduction 
electrons. In this paper we neglect this contribution because it is very small. The total 
susceptibility x is given by 
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d(M) 

X 



(14) 

As is well known, the Green's function Gfc CTO -(u;) is expressed by the retarded and the 
advanced Green's functions, G^ aa {u) and G{t aa {uj) as 

= ■ e(E* ka - fj) + GtM ■ 8(n - EL) (15) 

for \k\ ~ k F and u; ~ 0. [^(J E^ CT is the quasiparticle spectrum obtained by 

{^(4)} 4 = 0. (16) 
From the definition of above, we can show for |k| = k F that 



_d_ 

OH 



Ka = Ok(0) • E «kAf CT (0) [M5 M) M' + 7^7S kM M'(0) J «kA/v(0), (17) 
MM' V / 



where akMo-(^) is introduced by (§), and a k (a;) is the renormalization factor for the conduc- 
tion electrons, given by 



l/a k (cj) = "kMaM • \$MM' ~ j^kMM'M ) " "kM'a^)- ( 18 ) 

MM' \ dUJ ) 

Note that a k (to>) <C 1 for |u;| < T*. By use of (|i~5D and (|T7|), we can derive that 
9 / d \ 

— q— GkMM'(w) = (PkMM',M"M"'(u) ■ Sm"M"' ~ SkM"'M"(^) , (19) 

Co; M „ M „, \ Ccj / 

JTfzGkMM'{u) = ^ {</'kMM',M"M"'(w) + QkMM',M"M"'(w)} 
Cil M"M"' 

X ^Af "5jvf"M"' + ^S kM »/ M ''(^)J , (2°) 

where if and Q are defined by 

( PkMM',M"M"'( UJ ) = {GkMM"'(^)GkM"M'( UJ )} R ^ A \ (21) 

QkMM',M"M"'(^) = -2vria k (0) 2 5(u;)<5(£ k ) • a kA ./ (J (0)a k Af' (7 '(0)a kM // (T /(0)a kMW(T (0), (22) 

where the superscript R(A) represents retarded ( advanced ) function. On the other hand, 
we can show that 

lim lim G]<mm'"(^)G] 1 im"M'(m') = Q\lMM'M"M'"{u) + <P\iMM'M"M"'{u), (23) 

k— >k' lo^ui 1 

lim lim GkMM"'( w )Gk'M"M'(w) = (pkMM'M"M'"( UJ ), (24) 

ui^ui' k^k' 
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on the basis of the well-known FLT. p|,|20|,|2T[ By use of ([13]), (0), © and (H), we can 
prove that 

X = lim hrnXkH- (25) 



This relation is trivial whenever the magnetization operator is conserved P2"| , pTf but should 
be proved for our model. [0 Anyway, by use of (|25|), we can discuss our problem in terms 
of the more familiar FLT. (TT2|,^D,^T| So, \ can be expressed more explicitly as 



X 



Jf E / (^p( e ) + QM) kk M . (26) 

m + ^EJ^i {r* q fe <0 (<p q (0 + Q q (6')) m} 



M+— E p (e) 



where Tr represents the trace with respect to the angular momentum, and {M}miM 2 = 

Mi ■ $MiM 2 - Here, we have introduced the fc-limit and cu-limit of the four-point vertex 
r pp ',q(e, e';u>), and represents them as f pp ,(e, e') and f pp ,(e, e'), respectively : 

f pp ,(e,e') = hm hmf pp , q ( e ,e», (27) 

^ |q|— >0w— >0 

f- p ,(e,e') = hm hm f pp , q (e,e», (28) 

^ u>^0 |q|— >0 



where r pp 7 (e, e') 

ife(a;) 



PP' 1 

expresses the matrix whose (MM', M"M"') -component is V^, } MM , M /i M m(e, e'). They are 
expressed by Fig. 2. 

As is well known, r pp , (e, e') and r pp , (e, e') are related to each other by the following 
Bethe-Salpeter equation : [^0| . |T^1 

r^(e,e') = f- p ,( e ,e') + ^E / ^ {f pq (e, u>)Q^u,)i^{u>, e')} . (29) 

For instance, we can write the if-derivative and the cu-derivative of the selfenergy as follows 

-|;*pW = ^X:/^{^(«.«)W«)}. (30) 

J^pM = i £ / ^ {fj,( £ , W ) (4.M + 0„( u )) M) . (31) 
In this paper, we define the Van Vleck susceptibility xv an d the Pauli susceptibility xp 



as 
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Xv = lim Jim Xk(w), 

w-»Q |k|— >0 

Xp = X - Xv- 



(32) 
(33) 



According to d32[) and fl33|), we can divide x into xp and i n a unique way. The standard 
FLT tells us that xv — and x — Xp if t ne magnetization operator M is conserved, i.e., 
[H, M] = : [[H],|l2|] we briefly show this in Appendix A. Equations (|32] ) and (|33D can be 

expressed as 



Xv = ^E/^Tr{A^ p(e )A- (e )}, 
XP = ^E/^Tr{A-(e)g p (e)Aj( e )}, 

^ + ^E/^{fp q (^')^( e ')M} 

^ A P (e). ) . 



Ap(e) 



(34) 
(35) 
(36) 



M+ S-t 
OH 



-limit 



Equation (|34]) is depicted by Fig. 3. In the following sections, we investigate the expressions 



of x and XVi given by (|2q ) and fl34|) , further. 

Here we stress that the principle of division of the total susceptibility into (^) and 
([33]) is very clear ; the former does not contain any Q p (u) defined by (|22|) while the latter 
contains at least one Q p (uj) in itself. Equation fl35|) indicates that xp vanishes and the total 
susceptibility is given by xv m the case where Q p (e) = 0, i.e., where the Fermi surface 
disappears for some reasons or other. This is true for any multicomponent Fermi liquid 
systems. 



IV. MAGNETIC SUSCEPTIBILITY x° : IN CASE U = 

In this section, we calculate both Xv and x°p defined by (|34]) and ( |35"D in case U — 0, and 
show that our results are identical to those obtained previously by several authors. ||[7| 

X° = 4 E / ^ M [pI,mm,m'm>{u) + QI,mm,m'm>(u)}M', (37) 

N kMM' J 2711 1 J 
Xv = ~ E / ^%k,MM,M'M'(^)M' (38) 

a) Metallic case : For simplicity, we consider the case where the DOS for the conduction 
electrons is constant, i.e., p° c {oj) = p c . We also neglect the contribution from the spin of 
the conduction electrons, which is of order ~ p c /p ~ 0(V 2 /D 2 ) compared with the total 
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magnetic susceptibility. ( see §2. ). Moreover, we assume the six- fold degenerate /-electron 
spectrum and the spherical Fermi surface. The calculated results of (EJTj) and (1381) are 



X ° = 2g 2 J(J+l). V(0), (39) 
X v = 2g 2 J(J+l).Cp°(0), (40) 
A = 0.333, C = 0.152. 

Note that p°(u;), given by (|12"D, is the DOS for /-electrons. Moreover, the ratio x°p/Xv = 
(A — C)/C = 1.19 is derived. The Wilson ratio R° of this system, whose definition is given 
by (0), is given by 

R° = 1. (41) 

The obtained results ( |39"D ~ (f4~T|) are identical to those shown previously by several authors. 
@J^] Note that these results are not universal for J = 5/2 PAM in that it depends on the 
cu-dependence of p c {oj). 

b) Insulating Case : Next, we consider the insulating case where the Fermi energy p 
lies in the f-c hybridization gap. In this case, apparently Xp = because p°(0) = 0. So, the 
total magnetic susceptibility x° is equal to the Van Vleck susceptibility, which is given by 

xL = 2^J(J + l)-Cp (A_), (42) 
3V 2 

^ (A - )= >0 ' 

where A_ is the lower edge of the hybridization gap. In this case, x^ns remains finite and will 
take a large value. This result makes highly contrast to that of the insulating SU(6)-PAM, 
whose magnetic susceptibility is about V 2 /D 2 times smaller than fl4~2"|) . ( see Appendix A. ) 



V. MAGNETIC SUSCEPTIBILITY X IN THE D = oo LIMIT : IN CASE [7/0 

In this section, we analyze the magnetic susceptibility x i n the strong coupling regime, 
taking account of all the vertex corrections due to the strong Coulomb interaction U in a 
consistent way. In §5 and §6, we assume the six-fold degeneracy of /-electron spectrum and 
spherical Fermi surface in the extended zone scheme. Such a assumption will be allowed for 
our aim to elucidate the essential properties of the susceptibility under the influence of the 
strong Coulomb interaction. 

We study this problem concerning only the leading term with respect to 1/d-expansion ; 
we call such an approximation the d = oo approximation or the d = oo limit. In the d = oo 
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limit, the selfenergy becomes momentum-independent, as is well known. |17| For J = 5/2 
PAM, the electronic properties obtained by this approximation will be quite similar to that 



in the realistic d = 3 system. [18] So, the conclusions obtained by the d = oo limit analysis 
will be valid also for the d = 3 system. 



A. Brief review of the d = oo approximation 

In the d = oo approximation, the irreducible four-point vertex is considered as local, 
as well as the selfenergy. So, the reducible four-point vertex T is composed of the local 
irreducible four-point vertex T 1 and the non-local particle-hole Green's functions 0{uj), as is 
shown in Fig. 4. We explain the reasoning briefly : 

Fig. 5 shows the two examples of f/ 2 -order contributions for the susceptibility. We 
assume the d- dimensional hyper-cubic /-electron lattice for instance. Each site has z = 2d 
nearest neighbor sites, j, I, m and n represents the /-electron sites, and we have to take 
summation over them. Note that an /-electron on one site cannot hop to other sites without 
mixing with the conduction electrons. In Fig. 5 (a), (j,m) and (l,n) are connected by the 
Green's functions, while they are not in Fig. 5 (b). Here, we fix I and m as the £-th nearest 
neighbor sites. Such selections exist about ~ z^ ~ d^ different ways in total. In the spirit of 
1/d-expansion, an itinerant Green's function in the real space connecting the £-th nearest 
neighbor sites is scaled as ~ 0(d~^l 2 ). After the summation over j, n, I and m, while I and 
m are limited to the £-th nearest neighbor pairs, Fig. 5 (a) and (b) becomes of order ~ d^ 
and ~ d°, respectively. Thus in the d = oo limit, the irreducible four-point vertices, included 
in Fig 5 (a) for instance, are local and momentum- independent, as well as the selfenergy. 
On the other hand, the particle-hole pair should be regarded as non-local. Below, we take 
account of 0(d°) terms at most, i.e., the d = oo approximation. 

In this approximation, fc-summation becomes 

- Jpc°(tk)de k , (43) 

k 

where p c °(^k) is the unperturbed DOS for the conduction electrons. We see below that 
the results of this paper are sensitive to the value of p c °(0), but insensitive to the energy- 
dependence of p c °(efc) m heavy Fermion systems. So, we do not pay attention to the func- 
tional form of p c °(efc)- 
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B. The local Green's function, the selfenergy and the irreducible four point vertex in 

the d = oo case 



As is shown in the previous subsection, in the d = oo limit both the selfenergy and the 
irreducible four-point vertex are composed of the local Green's functions, g(u). When we 
assume the spherical Fermi surface, g(to) is given by 

J7 GkMM* M = g{^)' $MM> ■ (44) 

Thus, the local Green's function g(u>) is diagonal with respect to M, and it is independent of 
M in case H — 0. ( Notice that G^mm 1 ^) is proportional to the phase factor cx e^ M '~ M ^ fe , 
H which vanishes after the k-summation in case M 7^ M' . ) This is a remarkable simplifi- 
cation occurring in the d = 00 approximation. Thus, in the d = 00 limit, both the selfenergy 
and the four-point vertex become diagonal with respect to M even in a finite magnetic field 
H, and independent of M in case H = 0. 

For simplicity, we consider the case in the absence of ECF, i.e., E^ M = E$ for all M. In 
the absence of the magnetic field H, we can rewrite the /- electron Green's function given 
by ([|), into a simpler form as follows : 

GkMM'(u) = hkMM>Gk(u) + dkMM'G f (uj), (45) 

where 



-1 



-1 



(46) 



G k (u) = (co + fi - Ef - E(w) - 3V 2 /(u + p - e k ) 
G f (co) = (w + fi-E* 

hkMM' = Y, a V£_MoV'kM'o/'3V 2 , 
d-kMM' = VMM 1 ~ hkMM'i 

where E(u;) represents the selfenergy in the d = 00 limit, which depends on neither k nor 
M. 

Here, we consider the quasiparticle representation of G^{uj) given by (|46|) . For \uj\ < T*, 
the retarded Green's function is 

G?(") = Z t ] . f . ( 47 ) 
kV ' u)-E£ + i5 v ' 



( i.e., (-1/tt) ImG^(^) = Zl (u) ■ 5(cu - E$ 

where T* represents the characteristic energy, within which the quasiparticle can be well 
defined. ( T* corresponds to the renormalized Fermi energy. ) Here, Zi(u) is the momentum- 
independent renormalization factor of /-electrons, given by 
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, . , .j — ■ , I , (48) 

z{u)^{l-^i:^)) • (49) 

Both z(u) and Z\(oS) are the renormalization factors of /-electrons. Apparently, l/zi(0) = 
1/^(0) 3> 1 in the heavy Fermion system. 

Finally, in the d = oo limit, we can show that the irreducible parallel four-point vertex 
T I MM m , m ,(uj,uj') satisfies the following property : 

^MMM'M'i^j u') = T 7 (^, to') ■ 5 M M' } (50) 

i.e., V r MM M f M'{ u > ^0 i s diagonal with respect to M and independent of M ( see Fig. 6. ) 
This property is explained in Appendix. B, which is no more true for reducible vertices. 

Here, we summarise the results of this subsection : in our J = 5/2 PAM, in the d — oo 
limit, the angular momentum M is conserved in the local Green's function, the selfenergy 
and the irreducible four point vertex, because they are estimated as local processes in the 
d — oo limit. 



C. Vertex Corrections for the magnetic susceptibility 

Contrary to S(e) and T / (e, e'), the particle- hole Green's functions are never restricted to 
the local process, so (p(u) and Q(u) are not diagonal with respect to M. In d — oo, <p(u) is 
given by 

= 4Z^» ■ &iM = ^E^M ■ G k (w), (51) 

3 k 



which is a matrix with respect to M. In (plD, Gij(uj) ( G^lu) ) is the Green's function in 



the real- ( momentum- ) space representation. 

Here, we introduce a M- independent particle-hole Green's function in d = oo, ip a (u)), 

¥a{u) = — EE <PKMM,M'M'(v) ■ 1 
iV k M' 

= ( 1 -^ l )4E G k 2 (^) + ^ E G k(^)> (52) 

k k 

where 
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A = 1/3 = / ^ • frkMM- (53) 



<p a (w) is related with frequency-derivative of the selfenergy, which is depicted by Fig. 7 (a) 
Note that the summation on M' is taken in (|52]). 

Next, we also introduce another M-independent particle-hole Green's functions, <Pb{w), 



1V1 l iV k M' 



MM,M>M> (w) ■ M' 



+ 2 ^ + ^?(^)* W) + ^-^^(^T")^' (54) 

where 

B = 0.181 = / ^jr • Amm/M. (55) 
ikMM'= E ^kW^kM» CT MV k * M , v l/ k MvV(3^ 2 ) 2 . (56) 

ctct'M" 

<Pb(u;) is related with the magnetic field-derivative of the selfenergy, which is depicted by 
Fig. 7 (b). 

We also define the following functions, which are M-independent : 

Q a (u) = Yj E Qk,MM,M'M>(u) = -A ■ 2niz 1 (0)p(0)6(u) (57) 

iV k M' 

Qp(u) = | E E Qkmm,m'm> M • M'j = -5 • 27riz 1 (0)p(0)%) (58) 

Q 7 (w) = Q«(w) - Q^w) = -C • 27rLzi(0)p(0)«y(a;), (59) 
(i.e., C = A-B = 0.152), 

where we have used the relation dE^/de^ = a k (0) and zi(0) = a k (0) • 3V^ 2 /(/i — — S(0)) 2 
for |k| = fcp, which is rigorous because our selfenergy is k- independent. Note that p(0) 
is not an enhanced quantity. Especially, in the d — oo approximation, p(0) = p°(0) is 
satisfied rigorously on condition that E' is modified so as to fix the value of p, because of 
the momentum-independence of the selfenergy. 

Here, we notice that <Pa(^) 7^ ¥>&(<^) and Q a {ui) ^ Qp{u) in our model, which is the 
consequence of the fact that the magnetization M is not conserved, i.e., [H,M] ^ 0. This 
situation has been prevented us from treating this problem on the basis of the FLT. 

Then, we introduce the reducible four-point vertices r r (e, e') and r r, ^(e, e') so as to satisfy 
the following Bethe-Salpeter equations : 
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r (6, e') = r J (6, e') + f ^r 7 (e, u)cp r (u)T r (u, e'), (60) 

P*(e,e') = r'(e,e') + / ^V(e, u>)(<p r (u) + Q^u))^, e'), (61) 

where r = {a, a', b} and £ = {a, /?, 7}. For example, we can represent the mass-enhancement 
factor as follows : 

V 2 ) 



^(e + p-e,) 2 



= / ^r(c,e)^(e), (62) 
J 2m 

which is shown in Fig. 8. Here, we notice that J2 k fk,MM',M"M" — for M 7^ M' because 
of the phase factor, e^ M ~ M ^' Pk . In the same way, the Van Vleck susceptibility xv an d the 
total susceptibility x, given by fl34|) and (|26|) respectively, are expressed as 

Xv = 



-2/ J( J + 1) / |^ 6 (e) 1 1 + / |^r 6 (e, e>6 ( e0 1 , (63 ) 

X = -2^J(J+1) / il(^( e )+g / 3( e ))|l + / ^r^(e,e')Me') + g^(eO)}. (64) 

Notice that J2m'M" 1 1 ^ I Vmm,m i m'^ i ^Pm , m i ,m"M"M" = J2m> J J ^ I( Pmm,m'M'M' ■ T 1 ^ = 
M ■ J/rVbrVfe C M ■ JT b <p b and £ M M 2 = 2J(J + 1) for J = 5/2. In (H) ( (H and 
(0) ), we have omitted the contribution from the ( spins of the ) conduction electrons. We 
neglect it hereafter because its value is at most of order ~ p c (0)/p(0) ~ 0(V 2 /D 2 ) <C 1 



in the real heavy Fermion systems, [ 23 1 where p(0) ( p c (0) ) are the density of states of 
/-electrons ( conduction electrons ) at the Fermi energy, introduced by ( fl2|) ( by ([□]) ). 

In heavy Fermion systems, the behavior of Im<p a (e) or Im<pb(e) only for < — e <, T* 
plays a predominant role for the enhancement factor, as is discussed in Appendix C. Below, 
we show the relation between Im<p a (e) and Im<pb(e), which is derived in Appendix C. 

(i) In the metallic case : Here, we introduce the new particle-hole Green's functions, yv(e), 
as 

¥V(e) = <p b (e) - Q 7 (e). (65) 
Then, the following important relation is obtained : 

Im^(e) ^Im<p a (e) for < -e < T. (66) 

(ii) In the insulating case ( i.e., p(0) = ) : Here, we introduce vv( e ) as 
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<Pa'(e) = (f b (e) -Qy(e), 
g;(e) = -C ■ 2mz(A*_)p(A*_ ) ■ 5(e - A* ), 



(67) 
(68) 



where A* denotes the renormalized lower edge of the hybridization gap measured from the 
fermi energy, and p(A* ) = p c -(A* efc F ) 2 /3V 2 . Apparently, the relation p(A!_) = p°(A_) 
is satisfied very well. Because the relation l/z(A*_) = l/z(0) is naturally expected by the 
numerical calculation in Appendix D, the relation A* = z(0) ■ A_ is also expected. 
We can show that 

Im^(e) = lm<p a (e), for < -e - A*_ + < T*. (69) 

In the next subsection, we investigate the magnetic susceptibility \ and xv by use of the 
relation ( |66]) and (|69"D, respectively. 



D. The expression for the magnetic susceptibility in the d = oo case. 

In this subsection, we study the (|63"D and (|64D further and obtain the simple expressions 
for x, xv an d Xins- By use of the relations ^ = ip a ' + Q^ and (pb + Qp = W + Qa given in the 
previous section, we can check the following Bethe-Salpeter equations in a straightforward 
way : 

r 6 (e,e') = r a '(e, e ') + J ^'(e, a,)Q»r 6 (u;, e'), (70) 
T^(6,6') = r a '(e,e') + / ^T a \e,u)Q a (u)T^(u;,e'). (71) 

The first line is depicted in Fig. 9. Taking account of the relations (|66|) and f [59l) obtained 
in 55-3, we can estimate as 



r a (e,e') = T a \e,e>) for < {-e, -e'} < T*, 
T a > a (e, e') ^ r 6 >^(e, e') for < {-e, -e'} < T*, 

because the contributions to them from the quasiparticle excitations are dominant in heavy 
Fermion systems. These results play important roles in the following discussions. 
Then, we introduce the new enhancement factor, l/z'(e), as 

1 r df' 

^^^ss 1 *'*^ (73) 
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Considering the relation (|72|), 

l/z(e) l/z'(e) for < -e < < T* (74) 

is satisfied in heavy Fermion systems. In Appendix D, we show that the relation (f7^) is well 
supported by the numerical calculation by SOPT and SC-SOPT. 

Here, we introduce the enhancement factor for the magnetic susceptibility l/z H (e) as 

1 1+4,4^0 



z H {e) gMdH 



(75) 

11 H=0 

de' 



r fie 

I £ ( M(e - £/) + rM(e ' e ' )g « (e/) ) i^y (76) 



where (|76| ) is given for the metallic case. In the insulating case, T 6 '^, Q a in (|76| ) is replaced 
by T b , Q^, respectively. Apparently, 1/^(0) ~ ^/z(0) is expected. 

Then, we can rewrite the total susceptibility x> given by (0), and obtain a simple 
expression as 

X = -2g'J(J +1)1 1 ^-L- )Qa{e ) (2™S(e - £ ') + T^( £ , OQ.M) TT^y + 

= *°W?W +x5> ' (77) 

where x° — 2g 2 J(J + 1) • Ap(0), which is similar to the unperturbed value given by (p9|). 
And x'v i s gi ven by 



f°° de' , , , , , 1 r de d ^ , , . , 



-2 9 2 J(J+l)/^-( f )J-- (78) 
J 27Ti r(e) 

We can estimate that x' y = by use of (|74]), because 

f/f' L /' (/f (V 

k 

More accurate derivation of ( |T9"D is given in Appendix E. Equation (|77| ) ( and ( |64"D ) is 
depicted in Fig. 10. Thus, by use of ( |74|) and (£!|), we get the simple expression, 

This is the main result of this paper, and is proved to be rigorous when 1/z(cj) = Xjz'iuS) is 
satisfied identically. ( see Appendix D. ) 

Next, we also rewrite the Van Vleck susceptibility xv, given by fl63f) , as follows, 
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Xv = -2g 2 J(J + 1) J J ^ * Q 7 ( c ) (2vrW(e - e') + T% e')Q 7 (e>)) + X ' v , (81) 

where x'v * s gi ven by (ff8|). 

Before concluding this subsection, we consider the insulating case. In this case, appar- 
ently xp — and Xms = Xv at zero temperature. By use of the results obtained in the 
previous sections, we can show that 



X 



ins 



-2g*J( J + 1) / / ^-Lq< ( e ) (2-5(6 - e>) + F b ( e , (e')) ^ + Xv 
Xins ^(A*) z'(A*) Xy ' 1 J 



where x ms = 2g 2 J(J + 1) • Cp(Al), which is similar to the unperturbed value given by (|42|). 
is given by (^). Thus, we get the simple expression 

Xins — Xins ' ~ / a * \ • (83) 
Zh[A-) 

From the mathematical point of view, the obtained expressoin for x, given by (|77|) , is 
identical to (26). In fact, the energy-integration appearing in the definition of x-, given by 
(p6|), has been included in the enhancement factor at the Fermi energy, l/z#(0). Also, the 
expression (|8T| ) and (|82| ) are identical to (^) mathematically. However, the expressions 
obtained in this section show clear physical information : We can see that x and Xv ( 
or Xms ) are apparently enhanced by the interaction because they are proportional to the 
enhancement factor, 1/^(0) ( or l/z(Al) ). 



VI. FURTHER ANALYSIS FOR x , Xv AND xins 

A. Estimation for the Wilson ratio in the d = oo limit : in case U ^ 

In the previous section, we get the general expression for x and xv-, and find that they are 
enhanced by 1/z. Furthermore, we estimate the Wilson ratio in this section by examing the 
magnetic enhancement factor, 1/zh- For this purpose, we consider the T-linear coefficient 
of the specific heat and the charge susceptibility at first. 

(i) T-linear Coefficient of the Specific Heat ; 7 
In the metallic case, the 7 is given by, @ 

7 = 7°M(0)=7°M0), (84) 
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where 7 is the value in case U = 0. 7 is proportional to the mass-enhancement factor 
1/^(0) - l/z(0). 

(ii) Charge Susceptibility ; Xch 
By use of the relation 



9 v f ^ 



Q = ^ mH) (85) 



the charge susceptibility in the metallic case is given by 

+ / /^^^H 2 ^^ (86) 

where T a,a (e, e') is the antiparallel four-point ( reducible ) vertex. ( see Fig. B • 1.) The 
first term of ( |36"D turns out to vanish identically, as proved by (|79|) . 

Usually, in heavy Fermion systems, Xch is considerably suppressed by the strong Coulomb 



repulsion between /-electrons. [15] Here, we assume the following relation approximately : 



r a ' Q (0,0) = -T a ' a (0,0), (87) 

which is rigorous for J = 5/2 single-site Anderson model, and expected to be valid approxi- 
mately even for the lattice problem. ( In reality in the periodic system, the relation (^7]) may 
be modified by the magnetic correlation between different sites. But, we do not consider 
such an effect now. ) If we put Xch <C 1 in the strong coupling limit region, we get from 
(M) that 



-r a ' Q (0,0) • Az(0)p(0) = — = - for J = 5/2. (88) 

(iii) Magnetic Susceptibility ; x 
Here, we consider x of J = 5/2 PAM, in case there is no ECF and p(e) = constant for 
< -e < T*. At first, we consider the metallic case. By use of (|7§), (g) and r a ' a (0,0) = 
r 6 " 5 (0,0), the relation l/z#(0) = 1.2 ■ l/z(0) is obtained. Thus, x is given by 

X~1.2-X°-^. (89) 



So, the obtained Wilson ratio is R ~ 1.2 -R°. ( R is given by (M). ) In fact, the result 
is changed if the assumption ( j37|) is incorrect. Nonetheless, R ~ 1 is expected because the 
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contribution to R from the vertex corrections, T(0, 0) and T(0, 0), is about R—R° ~ 0.2 <C 1. 
p4| It should be stressed that the Wilson ratio for J = 5/2 PAM becomes equal to that for 
J = 5/2 impurity Anderson model, R = 1.2, if the magnetic correlation between different 
sites can be neglected. This result is nontrivial in the previous works, [|I|,|S|,|T2"fl and is 
contrastive to that for SU(2) PAM, R ~ 2. gS| 

On the other hand, because the relation < -r fe (0, 0) • Cz(0)p(0) < 1/5 is expected, the 
Van Vleck susceptibility flSlf ) is written by 

XV ~ ■ ^ (90) 



where = 2g 2 J( J + 1) ' Cp(0), which is similar to the unperturbed value, (pEPp. Thus, xv 
is enhanced also by 1/z. 

We stress that the RPA-type diagrams are included only in the polarization factor ( i.e., 
in R/R° ), and never included in the enhancement factor 1/z or 1/z' which brings the highly 
enhanced magnetic susceptibility observed in the heavy Fermion systems. ( see Appendix 
E. ) This situation also holds in the SU(N)-PAM. 

In conclusion, the Van Vleck susceptibility for the J = 5/2 PAM is strongly enhanced 
by 1/z, both in the metallic case and in the insulating case, and Xms for the J = 5/2 PAM 
is much larger than that in the insulating SU(N)-PAM without ECF. 



B. Consideration on the d = 3 system : in case U ^ 

Here, we briefly consider the magnetic susceptibility in the three dimensional system 
beyond the d — oo approximation. The relation x = X° ' (V z h){ z / V) is also derived in 
d = 3, because the equation ( |66D is also satisfied in d = 3 by the replacement of Q 7 (uj) with 
Q-y(w) ■ z(0)5(E%.) in d65|). On the contrary, the relation 1/z ~ 1/z' is not guaranteed in 
d = 3, because the relation ([50]) is rigorous only for the d = oo case. 

In the d = 3 system, there exist the correction from the non-local part of the irreducible 
four-point vertex, T^, which violates the relation fl5"0|). Fortunately, SOPT for J = 5/2 PAM 
by means of the 1 / <i-expansion shows that the correction term is approximately negligible 
in d = 3, so fl5T)|) will be satisfied very well. As a result, the relation Xv/Xv ~ Xp/x°p ~ V 2 
holds also in the three dimensional J = 5/2 PAM. 

At last, we make a more general but qualitative consideration on the original definition for 
Xv given by (|3~3]). Approximately, we can transform ( |3"1D to xv ~ ~ J-°oo dep(e) -^(A^ (e) z(e)) , 
where p(e) is the DOS for the /-electrons. ( see ( |ET| ) or (|E2"D. ) Here, we consider the quasi- 
particle contribution to Xv-, i- e - ; the contribution from the integration range, < — e <, T* . 
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It is not enhanced only when the relation |A w (e) • z(e)\ = 1 is satisfied. The enhancement 
of xv is n ot determined by the value of the enhancement factor at the Fermi energy, A w (0), 
but by the e-dependence of A aJ (e). ( From the definition of T*, \z(0)/z(— T*)\ <C 1 should 
be satisfied, and z(±oo) = 0. ) From the result of this paper, Xv/Xv ~ V- 2 * n d = 00 
J = 5/2 PAM, |A^(-T*)|/|A w (0)| < 0(1) should be realized. Such a property of A w (e) will 
not depend so much on the dimension of the system. Thus, xv should be enhanced by the 
strong correlation in three-dimensional systems. 



VII. DISCUSSIONS 

Here, we summarise the conclusions of this paper. At first, we obtain the general but 
abstract expression for the magnetic susceptibility on the basis of the orbitally degenerate 
FLT. In the next stage, we employ two simplifications : One of them is to assume no ECF 
and the spherical Fermi surface, and another is the d = oo approximation. Since our aim of 
this paper is to elucidate unambiguously the essential properties of the susceptibility under 
the influence of the strong Coulomb interaction, these over-simplifications will be allowed. 
Needless to say, it is significant to confirm the property for the zero ECF limit case. The 
opposite limit case, where only the lowest Kramers doublet contributes to the ground state, 



has already been discussed elsewhere. Pj2~8fl ( see Appendix A. ) 

After the two simplifications, we make further analysis on the expressions both for x and 
for xv, taking account of all the vertex corrections in a consistent way. Below, we summalize 
the results of this paper, which will hold qualitatively even in the three-dimensional case. 

(i) Metallic Case 

In this case the Fermi energy lies below the hybridization gap, which is a prototype of 
the ( Ce-compound ) heavy Fermion systems. [BJ The total susceptibility is given by x = 
limfc^oXfc(O) = Xv + Xp- In d = oo approximation, we get the simple expression, x = 
X°/zh{0). ( 1/zh(0) is the magnetic enhancement factor at the Fermi energy. ) We can 
also express x by use of the mass enhancement factor l/z(0) as x ~ 1-2 • X°/ Z (fyi which 
means that the Wilson ratio is R ~ 1.2. In the same way, the Van-Vleck susceptibility is 
also expressed as xv ~ 1-2 • Xv/ Z {®)- ^ n conclusion, both x and xv are proportional to the 
mass-enhancement factor. Our conclusion contradicts to the conclusion in ref. [|TJ. 

(ii) Insulating Case 

In this case the Fermi energy lies in the hybridization gap. This is a prototype of the so called 
Kondo insulators, some of which exhibit the large magnetic susceptibility experimentally at 
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T = 0. H In this case, xp — and the magnetic susceptibility Xins is given only by the 
Van Vleck susceptibility at zero temperature, we get the simple expression as Xins = 
Xins/zff^A^), where A* is the renormalized lower hybridization edge. Because 1/^(A* ) = 
1 / zh (0) is expected, the magnetic susceptibility for the orbitally degenerate model is strongly 
enhanced. This result is consistent with the result (i) because the Van Vleck susceptibility 
will be insensitive to the state of the Fermi surface. So, xv will be little affected by the 
superconducting transition. 

As our calculations are very lengthy and involved, we briefly summarise the mathematical 
analysis for the Van Vleck susceptibility, xv '■ As is shown by (^), Xv is given by the cu-limit 
of the dynamical magnetic susceptibility. After the d = oo approximation, it is written as 
Xv ~ + rVf>); where ifb and Y b are given by (^) and (|60|) , respectively. ( Here, the 
symbols for the energy-integration are implicit. ) On the other hand, l/z(u) ~ (1 + T a ip a ). 
T a and T b are related by the Bethe-Salpeter equation given by (|70|) , T b = T a + T a Q 1 T b . ( 
Here, we have identified r a with T a ' for simplicity. ) By use of this equation, we can finally 
show that xv ~ / duj{l/z(uS) ■ Q 7 (uj) ■ l/z(u)} ~ l/z(0). In the same way, x ~ 1/^(0) i s 
derived. 

Here, we make the physical consideration on xp an d Xv i n terms of the one-body 
picture, which was done previously by Anderson and Zou. ||TT| In the one-body picture, 
Xp/Xp ~ 1/2(0) an< l Xv/Xv ~ (^ ~ f 1 )/ {E* ~ lAi where z(0) is the renormalization factor 
at the Fermi energy and E$ is the renormalized /-electron spectrum. In the mean-field 
approximation, where the frequency-independent renormalization factor z cons t is assumed, 
then is strongly renormalized towards the Fermi energy so that Xv/Xv ~ l/^const- On 
the other hand, from the viewpoint of the FLT, renormalization is caused by the strong 
energy-dependence of the selfenergy. The renormalized value of critically depends on 
the energy range of the coherent region ( around the Fermi energy ), in where the energy 
dependence of the selfenergy is large. Thus, the question Ho what extent xv is enhanced 
or not by the strong correlation ' is never trivial within the one-body picture. The simple 
mean-field approximation never answer this question. 

Our analysis shows that the excitation of the quasiparticles, which are well defined only 
within T* , bring the enhancement of xv, Xv/Xv ~ 1/^(0). Our work also suggests that 
E-f is strongly renormalized to the Fermi energy. 

Finally, we point out some future problems. At first, it is interesting to estimate the 
influence of the electronic crystal field, or the shape of the Fermi surface ( i.e., the shape of the 
lattice ). These effects are ignored in this paper. Secondly, the effect of the antiferromagnetic 
fluctuation should be taken into account correctly. The antiferromagnetic fluctuations will 
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make the magnetic susceptibility x smaller than our prediction, although they cannot change 
the value of x an d Xv drastically. In fact, some heavy Fermion compounds are under the 
influence of the prominent antiferromagnetic fluctuations, which may be the driving force 
to the superconducting state. Such a study will give us much information on the electronic 
properties of the heavy Fermion systems. 
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APPENDIX A: THEORY OF THE MAGNETIC SUSCEPTIBILITY FOR 

SU(N)-PAM 

In this Appendix, we briefly consider the magnetic susceptibility for SU(N)-PAM. In case 
U = 0, it is easy to show that Xv * s mucn smaller than x°p- Here, we study the case U ^ 0, 
making no use of the d = oo approximation. The Hamiltonian without ECF is given by 

#o = E + E E f fLhM + J2( v fL^M + voUhM), (ai) 

kM kM Mk 

Hi = — E Zk-qM/k'+qM'/k'M'/kM- (A2) 
Z kk'qM^M' 

We note that the conduction electrons have the six-fold degeneracy. This model is analyzed 
usually by the slave boson technique. [p6f 



Here, we assume the existence of the magnetic field H along 2- axis. At first, we assume 
that both E? and appearing in (|A1 ) are shifted by the same Zeeman energy as 



E f M = Ef+g ff x B M-H, 

(Ad J 

ekA/ = £k + 9c^bM ■ H. 
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Only in the special case Qf — g c , the magnetization operator M is conserved and \ is derived 
after the Luttinger's manner. [13|;[15| But, we consider the case g/ ^ g c here. The Green's 
functions for /-electrons G^mO^) and that for the conduction electrons G% M (u) are given by 



G kM (u) = ( u + ii - Ef - E kM H - ^— - ) , ( A4) 

GUM = G kM (") ■ V 2 /(u e k ) 2 , (A5) 

where [uS) is the selfenergy. 

Here, we define the three types of the enhancement factors. 

' 1 - ^-£ fcM M, (A6) 



' l + ^W. (A7) 



z H (u) g f MdH 
1 1 



z'{u) z H (u) 



(A8) 



-limit 



where (|A8|) is defined by (|36j). From the definition of Xp, given by (|35|) , we get 

f o (i g c v 2 \ ( i v 2 r 1 f i ^ . rAq , 

Xp Xp s/0"-OV Uo) (a*-c*,)V Wo) ff/ (a* - c*,) 2 1 

M 

where p(0) is the DOS of the /-electrons at the Fermi energy. On the other hand, by use of 
(0), we get 

1 2 , ,2 f°° ( 1 9c V 2 \ . ( 1 V 2 \ 



g f (uj + fi - e k ) 2 J \z(u) (u + fi-e k ) 



2 



x^GW(aO + |^(w)V (A10) 

In case = <7 C , l/z(u) = l/z'{u) is satisfied rigorously. Then, xp = Xp/2h(0) and 
Xv = is derived. Although l/z(u) = 1/z'{uj) is not satisfied rigorously in case gf ^ g c , 
its discrepancy will be of order ~ (V/D) 2 /z at most because the relation a ^ H ^ £m(0) ~ 
^E M (0) • D 2 /\/ 2 is expected from (§2]) and ref. 

Thus, xv is negligible compared with \p- So, the susceptibility for the insulating case, 
Xins, cannot become large. This reasoning has been already applied to the system with the 
Kramers doublet ground state in the strong ECF limit. [p|,p8[ In the Gutzwiller approxi- 
mation, Xv/Xv i n this model is not enhanced because it is proportional to the conduction 
electron DOS, which is not renormalized. But the more detailed analysis should be re- 
quired. Finally, we point out that xv — for any model where M is conserved, because 
1/z{uj) = \/zh{oj) is satisfied then. 
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APPENDIX B: SOME PROPERTIES OF THE IRREDUCIBLE VERTEX 



At first, we show that we can ignore the Pauli's principle in (|j). The term excluded in 
(|j) due to the Pauli's principle is, in the real-space representation, given by 

U fLlfiMfmfiM, (Bl) 

{i},M 

where {i} represents the set of the /-electron sites. Apparently, (|B1|) represents a constant 
energy shift and is independent of M in the paramagnetic state. Therefore, we can ignore the 
Pauli's principle in the following discussion if all the diagrams are taken into consideration. 
In other words, contributions from the diagrams violating the Pauli's principle cancel out in 



each order of the perturbation in the paramagnetic state. [27 

Then, we investigate the general properties of four-point vertices. We classify four-point 
vertices into the parallel vertex T and the anti-parallel vertex T, as is shown by Fig. B 
• 1. Needless to say, T never contributes to -J^Y< m (uj). Moreover, ^T, M (u) also have no 
contribution from T because Pauli's principle is ignored now. 

In d = oo limit, by neglecting the Pauli's principle, we can regard that the irreducible 
parallel four-point vertex T 1 satisfy the relation 

r MM,M'Af( e ) e ') = r/ (e, e) ■ 5 mm', (B2) 

because the local Green's function is independent of M and diagonal with respect to M. 

On the other hand, when we investigate -^-Hm(^)i we have to take account of the con- 
tribution from both T and V. The irreducible antiparallel four-point vertex T 1 satisfy the 
relation 

T I MMM'M'^t') = T\e,t') for any M. M'. (B3) 
by neglecting the Pauli's principle. 



APPENDIX C: THE CONSIDERATION ON <pa(B)(«) and the PROOF OF (66) 

AND (69) 

At first, we consider the following calculation as a preparation of calculating the en- 
hancement factor. 

00 dtii 

—r'(0,u;)<p a{b) (u)r'(u;,0), (CI) 

-oo ZTT1 
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where T'(e, e') is an irreducible four-point vertex with respect to the Coulomb interaction U. 
Then, 

{eq. (ci)} = ^r'(o )W )^Hr'( Wl o) 

= jf° ^ {lm{y, a(fe) (C)}Re{r' 2 (C,0)} + Re{^ (fe) (C)}Im{r' 2 (C, 0)}} . (C2) 



In the first line of the r.h.s. of (C2), we have changed the path of w-integration as shown in 
Fig. CT. Even if T f is replaced with reducible vertex T, this procedure is correct because ip a 
~ (pf, ~ cij~ 2 for uj ^ D. The last term in ( |C2|) can be obtained by substituting the following 

»(6)( 



spectrum representations ( which is possible because <p a (&\(±oo) = r /2 (±oo,0) = 0.) 



= - / ^ {^^- + ^^-} ' I'^VO. (C3) 
J ir [w-( + if| uj — Q — it] j v; 

I>,0) = - ! d M^P- + ■ Imr« 2 (C,0), (C4) 

and by performing ^-integration at first. In the last line of (|U2|) , only Re{r'(£,0) 2 } for 
< — ( < T* contributes predominantly to the value of the (^-integration, which is largely 
enhanced. Thus, in heavy Fermion systems, the behavior of lm(p a ^) (e) only for < — e <, T* 
is dominant in determining the enhancement factor. 

In the next stage, we prove the important relation given by (|56]). The spectral represen- 
tations for the local Green's functions, given by flip) , are 



3 ( 'V)=/dc{^L + ^U- /><"«), (C5) 

J [UJ — C + 1? ? LU — Q — 17] ) 

P(0 = --Im^£G* MM (C) = --lm^(C), 

7T iV k 7T 

P / (C) = --Im^E^MMG'f(C), 
n N k 

where p^(C) represents the DOS of the localized /-electrons, and dkMM is given by (|46|). 
Apparently, > p^(C) for any ( and p^(0) = 0. On the other hand, we can also express 
the particle-hole Green's functions in the spectral representation as follows ( see ( |C5| ) ) : 

where in deriving the last line above, we have done the partial integration. 
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At first, we investigate the case where the system is metallic at T = 0, i.e., p(0) 7^ 0. ( 
see Fig. 1(a).) For < -u < T*, 

1 W 2 

k (oj + p — E? — h(tj)) 

W 2 

= PcM — o. (C7) 

m ; (w + p - Ef - E(w)) 2 1 ; 

So, we obtain 

9 , s / \ 1 Id,. 3V 2 
—piuj) ^ -2p(u) —. — ■ + -piw ? . C8 



Considering that the first term of the r.h.s. of (|Cq ) has the enhancement factor 1/z contrary 
to the last term, we can neglect the last term of (08) for < —uj <, T* in metallic heavy 
Fermion systems. On the other hand, pf (uj) is non-zero only for uj ~ \E?* — p| > 0, 
corresponding to the renormalized Ef spectrum, and uj ~ ±U /2, (U T*), corresponding 
to the broad satellite on both sides of the Fermi energy H. Especially, p / (0) = 0. Thus, 



both \mG k (uj) and ImG{ (uj) make little contribution to the behavior of (|66D or (|69"1), so 
they are negligible. By use of the relation for < — uj < T*, 

'4 £(-i/»ggM ■ g , + „ -Ef - w ■ < C9 > 



we can show from (p4|) that 

Im^H S Im |sl £ G k 2 (^) + (A-B)(^Y, <%(<*) ~ 2niz(0)p(0)5(uj) 
= lm{<p a (uj) + Q 7 (u)}, 



for < — uj < T* . Here we have used (|U6|) and (|C8|). We stress that when U = and 
p°(u;) = constant with respect to uj, the relation fl66|) is rigorous for uj < 0. 

In the same way, we consider the case where the system is insulating, i.e., p(0) = 0. ( see 
Fig. 1(b). ) We prove the important relation given by (j69|) . Note that p(uj) = for A* < 
uj < Ef* — p, where A* and Ef* are the renormalized lower edge of the hybridization gap and 
the renormalized local /-electron level, respectively. They are given by A!_ = A_ ■ z(0) and 
Ef* ^ (E f - p) ■ z(0) + p. O In this insulating case, taking account of the fact lmE(u) = 



for A* < uj < Ef* — p, we can show that 



V/ 2 

pM = PM— T( ■ 9(A*_ - uj), (C10) 

(uJ + p — Ef — 11{U))Y 



for A*_ — T*<uj< Ef* — p. So, we get for uj < 0, 
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9 ( \ 



{eq.(C8)}-p(A*)5(u;-AI). 



(Cll) 



Thus, for < — uj < — A*_ + T*, we can show the relation ( |69|) in the insulating case. 
Here, we comment that if the anisotropy of the Brillouin zone is taken into account in d = 3 
system, the step function in ( |C10| ) becomes a continuous function because of the van Hove 
singularity. So, the delta function in ( |C11D has finite ( but narrow ) width. In this sense, 
the relation (|69| ) is less universal than the relation (|66|). 



APPENDIX D: NUMERICAL CALCULATIONS FOR 1/Z(u) AND l/Z'(u) 

In this Appendix, we explain the method and the results of the numerical calculation 
by SOPT and SC-SOPT with respect to U in the d = oo limit. We calculate both l/z(u) 



and l/z'(u) and check the relation ([74]), which is derived in the analytical way. The relation 
([74]) is very significant because our main results (|80D and (]83|) are based on it. 

At first, we calculate l/z(u) and l/z'(u) by SOPT, which are depicted in Fig. D-l(a) 
and D-l(b), respectively. Here, we use the constant DOS for the conduction electrons, i.e., 
P° c {oj) = p Q c . The numerical results of their real part are shown in Fig. D-2. The obtained T* 
in SOPT is about \E? — p\. Surprisingly, the resultant l/z(u) and l/z'(u) are almost the 
same not only for \uj\ < T* but also for ^> T*. But, this result is modified when p° c {oS) 
has drastic energy dependence for < — uj < T*. By use of SOPT, we also calculate l/z(u) 
and l/z'(u>) for insulating case. We put Fermi energy p in the middle of A_ and E* . The 
numerical results of their real parts are shown in Fig. D-3. In the framework of the SOPT, 
A_ is renormalized to be A* = z(0) ■ A_. 

In the second stage, we calculate l/z(u) and l/z'(u) by SC-SOPT. The irreducible 
four-point vertex T 1 is composed of at most two local Green's functions, and is shown by 
Fig. D-4. In SC-SOPT, we construct l/z(u) and \jz\uj) with T 1 and g{oS) by referring to 
their definitions, flB2p and ([73]), respectively. ( Then, [/-linear terms of V 1 turn out to give 
no contribution, as is shown below. ) The numerical results of their real parts are shown 
in Fig. D-5. Here, we assume p° c (uj) oc (u + D) 2 . The obtained T* is renormalized to be 
~«(0) • \Ef + E(0) - p\. 

Surprisingly, in spite of p° c {oj) ^ p° c) the resultant l/z(u) and 1/ z{u) are almost the same 
not only for |u;| <, T* but also for |u;| ^> T*. The discrepancy becomes almost invisible 
in case p° c {uj) — It becomes large only when p° c {ui) has drastic energy- dependence for 
< — uj <T*. In real heavy Fermion systems, T* is renormalized as is pointed out by 
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SC-SOPT, so the susceptibility will be quite insensitive to the energy dependence of the 

In calculating for \jz\uS) in case p° c {uj) ^ p° c) we do not fix the value of C to ([59]), but 
adjust it so as to the equation ([T^) is satisfied. Then, the obtained l/z'{u) has no RPA-type 
enhancement, as well as l/z{oS). The resultant values of C ca i are given by Table D • I. The 
larger the l/z(0) is, the smaller \C — C ca \\ becomes, as is expected in Appendix C. ( see (|C5|). 
) Taking account of the smallness of 1/z in our calculation, we can say that our numerical 
result supports (|7"4|) quite well even in case p° c {u) ^ p° c . 



u 


Ccai 


l/z(0) - 1 


0.5 


0.143 


4.67 


0.3 


0.140 


3.64 



Table D-I : The values of C ca i which satisfies by eq. (|79|) by SC-SOPT. All the parameters 
except for U are the same as that used in Fig. D-5. 

Here, we do not calculate the magnetic enhancement factor l/z H (u) = 1 + ^t^£ m (u;) 
because l/z H (u) calculated by SC-SOPT will contain much contribution from RPA-type 
diagrams. We know that their contribution is limited only to the polarization factor, Rj R° ~ 
1.2. ( see §6 and Appendix E. ) On the calculation of 1/ zh{uo) by SC-SOPT, we have to take 
caution against the unphysical enhancement of the Wilson ratio caused by the RPA-type 
diagrams. In the rigorous perturbation calculation, most of them turn out to be canceled 
with other diagrams not contained in SC-SOPT, and amount to a little contributions totally. 

APPENDIX E: PROOF OF (79) 

In this appendix, we examine the relation \v = 0, shown by ([79]), in detail. To do this, 
we study the following integrals at first : 

1 E r £-Gl(e)/z'(e) = 1 £(-2i)Im f ^G k (e) . §-{ Zl (e)/z'(e)}, (El) 

1\ k J-oo 2,711 iV k J-oo Z7T1 Oe 

l E T ^GC(e)/z'(e) = 1 E(-2i)Im f £-M(e) ■ §-{z(e)/z'(e)}, (E2) 
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where we have done partial integration. Because J^(z(i)(e)/V(e)) <C 1 for < — e <, T* is 
expected, both ( [El]) and (|E2|) are not enhanced. As x'v * s expressed by both ( [El] ) and ( |E2| ) 
( see ( [66]) and (|78|) ), we conclude that Xv * s no ^ enhanced, and small quantity. Even if 
Xv 7^ 0, a slight modification of the value of C in the definition of Q^(uj) can make x'v = 0- 
( As is shown in Appendix D, this is confirmed well within SC-SOPT. ) This reason is 
as follows : the behavior of z'(e) for < — e <, T* is sensitive to the value of C, and the 
absolute value of ( |E1| ) is largely enhanced when ([F4|) is not satisfied. Thus, the relation (|79|) 
is justified. 

1/ ' z'{uj) can be rewritten as 

l/z'H = (l +J ^Ti a .U^e) (Paf (e)\ ■ (l + Ux' v (3/2gij(J+ 1))) , (E3) 

where r^ n _ tad (d;, e) is derived from T a '(tu,e) by dropping the reducible terms with respect 
to U . When we choose the value of C so that Xv = is satisfied, ( which will be close to 
the value given by (|59]), ) the enhancement factor l/^'(0) has no contributions from any 
irreducible vertices with respect to U, i.e., from any tadpole diagrams. ( see Fig. E-l. ) 
Moreover, the mass-enhancement factor l/z(0) also has nothing to do with any tadpole 
diagrams, that is, 

1AH = 1 + / |^ on _ ta >, eK(e). (E4) 

Then, we can conclude that both l/z(u) and l/z'(u) has no contributions from RPA-type 
enhancement. This results suggests that the relation l/z(u) = Xjz'iiS) holds for a wider 
range of uj. 
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Figure 

• Fig. 1 : Effective band structure for J = 5/2 PAM. E? is the four-fold degenerate 
local /-electron spectrum, and E£ is the two-fold degenerate quasiparticle spectrum, 
respectively, (a) Fermi energy /j, lies below A_ + /i. The system is metallic, (b) /x lies 
between A_ + /i and . The system is insulating. 

• Fig. 2 : (a) Four point vertex rJ%- M " M "'(e, e'; u). (b) fc-limit ( or cu-limit ) four 
point vertex r£ ,(e, e'). 

• Fig. 3 : The schematic structure for xv- 

• Fig. 4 : The structure of the vertex correction for the ( effective mass or magnetic 
) enhancement factor in the d = oo approximation. Here, I represents the local 
irreducible four-point vertex, T 1 , and <p represents the particle- hole Green's functions, 
connecting between two sites. /, m, n and j represent the /-electron sites, on which 
we have to take summation. 

• Fig. 5 : Two examples of £/ 2 -order contributions for the susceptibility, j, I, m and 
n represent the /-electron sites, on which we have to take summation. 

• Fig. 6 : ^mm M'M'{ e -> e ') represents the irreducible parallel four-point vertex with 
respect to the particle- hole pair. In the d — oo limit, it is diagonal with respect to M 
and independent of M. 

• Fig. 7 : The definition for two kinds of the particle-hole Green's functions, tp a (uo) 
and tpb(u). Both of them are independent of M. I represents the /-electron sites, on 
which we have to take summation. In (b), factor M' comes from the Zeeman term. 

• Fig. 8 : The schematic structure for the energy- derivative of the selfenergy. This 
is composed of { T 1 , (p a }. 
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Fig. 9 : The first line represents the Bethe-Salpeter equation relating between T a ' 
and T b . The second line is derived by using the first line twice. 

Fig. 10 : The schematic structure for xv- This is composed of { T 1 , ip a i }. We 
stress that two enhancement factors appear in total on both sides of Q~,(ui) in the last 
line. Xv turns out to vanish identically. 

Fig. B-l : r(e, e') represents the parallel four-point vertex and T(e, e') represents 
the anti-parallel four-point vertex, respectively. 

Fig. C-l : Complex integration path C in eq. (|C2|). 

Fig. D-l : The diagrams for the enhancement factors in SOPT. The broken lines 
represent the Coulomb potential U. (a) {l/z(u) — 1} in SOPT. This is constructed 
by the unperturbed (p° a {oj). (b) {l/z'(u) — 1} in SOPT. This is constructed by the 
unperturbed <£>°,(co>). 

Fig. D-2 : Numerical results for the frequency dependence of the real parts of 
enhancement factors by SOPT, for the metallic case. The line and the broken line 
represent {l/z(u) — 1} and {l/z'(u) — 1}, respectively. The energy region for large 
enhancement in SOPT is about <\fi-Ef\. We put U 2 = 1, 3V 2 = 0.16, E f = -0.3, 
kp = 0.8857T and e k = 2(|&;|/7r) 3 — 1, respectively. 

Fig. D-3 : Numerical results for the frequency dependence of the real parts of 
enhancement factors by SOPT, for the insulating case. The line and the broken line 
represent {l/z(u) — 1} and {l/z'(u) — 1}, respectively. We put U 2 = 1, 3V 2 = 0.25, 
E f = -0.5, fi = (A_ + Ef)/2 and e k = 2(|/c|/7r) 3 - 1, respectively. 

Fig. D-4 : The diagrams for the irreducible parallel four point vertex r 7 (e, e') used 
in SC-SOPT. 

Fig. D-5 : Numerical results for the frequency dependence of the real parts of the 
enhancement factors by SC-SOPT, for the metallic case. The line and the broken line 
represent {l/z(u) — 1} and {l/z'{u) — 1}, respectively. The energy region for large 
enhancement in SC-SOPT is about ,< z(0) ■ \fi — E* |, which is renormalized compared 
with that by SOPT. ( see Fig. D3. ) We put U 2 = 0.5, W 2 = 0.16, E f = -0.3, 
kp = 0.8857T and e k = 2(\k\/ir) — 1, respectively. 

Fig. E-l : The diagram for 1/z' (u). The first term is irreducible and the second term 
is reducible with respect to U. As the second term includes Xvi ^ vanishes identically. 
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